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We present an analysis of the semi-inclusive decays Bs —>■ Df Xi^u and Bs —>■ D'l~XI^v, where 
X denotes a final state that may consist of additional hadrons or photons and I is an electron 
or muon. The studied Bs decays are contained in the 121.4 fb“^ T(5S) data sample collected 
by the Belle detector at the KEKB asymmetric-energy e"'"e“ collider. The branching fractions of 
the decays are measured to be B{Bs —>■ Df Xt^v) = [8.2 ± 0.2(stat) ± 0.6(syst) ± 1.4(ext)]% and 
B{Bs —>■ Ds~X£^u) = [5.4 ± 0.4(stat) ± 0.4(syst) ± 0.9(ext)]%, where the first two uncertainties are 
statistical and systematic and the last is due to external parameters. The measurement also provides 
an estimate of the Bi*'^ production cross-section, a{e^e~ —>■ Bi*'^= (53.8 ± 1.4(stat) ± 
4.0(syst) ± 3.4(ext)) pb, at the center-of-mass energy ^/s — 10.86 GeV. 

PACS numbers: 14.40.Nd, 13.20.He 


INTRODUCTION 

Analyses of semileptonic decays B —>■ X^tv, where 
X(. denotes a hadronic final state with a charm quark, 
play an important role in the determination of the CKM 
matrix element \Vcb\- The extraction of jUcbl from the 
measured decay rates relies on form factors that describe 


the accompanying strong interaction processes. Measure¬ 
ments of semileptonic Bg decays provide complementary 
information to test and validate the QCD calculations 
of these form factors. Since large Bg samples have be¬ 
come available at Belle and the experiments at the Large 
Hadron Collider, the interest in the topic of semileptonic 
Bg decays has intensified recently. Theoretical predic- 



3 


tions of form factors and branching fractions are based on 
QCD sum rules [1-4], lattice QCD [o, 6] and constituent 
quark models [7-15]. The predicted exclusive branching 
fractions vary from 1.0% to 3.2% for Bg —^ Dgtv decays 
and from 4.3% to 7.6% for Bg —>■ D*giv decays. There are 
also predictions for the modes with higher excitations of 
the Dg meson, denoted hereinafter by “£)**”. The LHCb 
and D0 experiments have measured the semi-inclusive 
branching fractions of the decays Bg —)■ Zlsi(2536)X/x+i/ 
and Bg —>■ Zl*2(2573)X/r+z/, where the D** mesons were 
reconstructed in D^*^K final states [16, 17]. The inclu¬ 
sive semileptonic branching fraction of Bg —>■ Xiv de¬ 
cays was recently measured by Belle and BaBar [18, 19] 
and found to be in agreement with the expectations from 
SU(3) flavor symmetry [20, 21]. We report here the first 
measurements of the semi-inclusive branching fractions 
B{Bg —>• DgXiv) and B{Bg —>• D^Xiv) using the Belle 
T(5S') dataset. The number of B'f''Bi*'^ pairs in the 
dataset, 

Nb^b. = ^(e+e- ^ sWsi*)) • Cri^s ), (1) 

where (T{e'^e~ —> Bi*'^Bi*'^) is the production cross- 
section and Cr{5S) is the integrated luminosity, is the 
limiting systematic uncertainty in this measurement and 
other untagged Bg measurements at Belle [22] . The value 
Ng^Bs — Cf.lil.S) X 10® was obtained from a measure¬ 
ment of the inclusive Dg yield in the dataset [23]. The 
measured Bg —>■ DgXiv yield, together with an estimate 
for the branching fraction B{Bg —>■ DgXli/), provides an 
alternative way to determine Nb^Bs- ^ similar approach 
was already pursued by the LEP experiments [24-26] and 
LHCb [27]. 

DETECTOR, DATA SAMPLE AND SIMULATION 

The Belle detector located at the KEKB asymmetric- 
energy e+e“ collider [28] is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector, 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel¬ 
like arrangement of time-of-flight scintillation counters 
(TOE), and an electromagnetic calorimeter (ECL) com¬ 
prised of CsI(Tl) crystals located inside a superconduct¬ 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside of the coil is instrumented 
to detect mesons and to identify muons (KLM). The 
detector is described in detail elsewhere [29]. 

This analysis uses a dataset with an integrated lu¬ 
minosity of jCt(5S) = (121.4 ± 0.8) fb“^ collected at a 
center-of-mass (CM) energy of ^/s = 10.86 GeV [30], cor¬ 
responding to the mass of the T(55') resonance. The Bg 
mesons are produced in pairs in the following production 
modes, with the respective production fractions given in 
parentheses: B*gB* ((87.8± 1.5)%), B*gBg ((6.7±1.2)%) 


TABLE I: Nominal branching fractions of D** decays to dif¬ 
ferent final states in the MC simulation. The branching frac¬ 
tion of the Dg —>■ DgX decays is set to 100% and this crossfeed 
is included in the calculation of the branching fractions to the 
DgX final state. 


Y Y - 

Branching fraction [%] 

DsX Y D*X Y 

£,*0(2317) 

100 

63 

0 

£,i(2460) 

100 

3 

0 

£,i(2536) 

0 

0 

100 

£,*2(2573) 

0 

0 

100 


and BgBg ((2.6 ± 2.6)%) [31]. All production modes 
are considered for the analysis. Moreover, we use a 
62.8 fb“^ sample collected below the production thresh¬ 
old for open B production to study the continuum pro¬ 
cesses e+e” qq (q = u, d, s, c). 

A sample of simulated events with a size correspond¬ 
ing to six times the integrated data luminosity is gener¬ 
ated using Monte Carlo (MC) techniques. The simulated 
data emulate the different types of events produced at 
the T(5iS') CM energy, comprising events with B and Bg 
decays, bottomonium production and the qq continuum 
processes. The events are generated with the EvtGen 
package [32] and are processed through a GEANT [33] 
based detector simulation. Einal state photon radiation 
is added with the PHOTOS package [34]. 

The branching fractions in the simulation are set to 
the latest averages from the Particle Data Group [31]. 
However, for semileptonic Bg decays only measure¬ 
ments of the Dgi{2536) and 11*2(2573) modes are avail¬ 
able, so we use instead the calculations from Faustov 
and Galkin [8], who predict the full set of branch¬ 
ing fractions and thus provide a self-consistent picture 
of the semileptonic width. The Bg semileptonic de¬ 
cay modes considered in this analysis, with their corre¬ 
sponding branching fractions given in parentheses, are: 
Dg (2.1%), D*g (5.3%), £1,1(2536) (0.84%), £*0(2317) 
(0.36%), £,i(2460) (0.19%) and£;2(2573) (0.67%). The 
decays Bg —>• D'f'^f.v are described by the Caprini- 
Lellouch-Neubert model [35], based on heavy quark ef¬ 
fective theory [36]. Assuming SU(3) flavor symmetry, 
the form factors of the semileptonic Bg decays are taken 
to be identical to the ones measured in the correspond¬ 
ing B decays [37] ; we use the following values of the form 
factor parameters: po = 1.186 for Bg —>• Dgiv decays, 
and p = 1.207, Ri = 1.403, i ?2 = 0.854 for Bg D^iv 
decays. The —>■ Dg*£v decays are described by the 

Leibovich-Ligeti-Stewart-Wise (LLSW) model [38] orig¬ 
inally developed for B —>■ D**£v decays. We replace in 
this model the B and £** masses by the Bg and £** 
masses, respectively. The nominal branching fractions 
for the £** decays in this analysis are listed in Table I. 
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ANALYSIS OVERVIEW 


This analysis is based on samples of reconstructed 
Dji~^ and pairs [39]. Incorrectly reconstructed 

Dg and D* candidates constitute a large background in 
the analysis. We therefore perform fits to the Ds mass 
distributions to determine the yields of events with cor- 
rectly reconstructed ' mesons. These events contain 
the following signal and background categories: 

1. e“''e“ —>■ cc continuum; 

2. B ^ Ds K£v decays, which have a branching frac¬ 
tion of (6.1 ± 1.0) X 10-4 [31, 40, 41]; 

3. opposite-B(^s) background, where a lepton candi- 
date is combined with a Ds meson from the sec¬ 
ond Bg in the event; the lepton candidate can be 
either a primary lepton from a B(^s) Xtv decay, 
a lepton originating from a secondary decay or a 
misidentified hadron track; 


4. same-B(^s) background from secondary leptons and 
from hadron tracks misidentified as leptons, which 
stem from the decay of the same Bg meson as the 
reconstructed Ds meson; 


5. signal: in the DgXlv channel, the signal comprises 
Bg Ds£y decays and crossfeed from Bg 
and Bg —>• Dl*(.v decays; in the D*X£iy channel, 
the dominant signal contributions are Bg —>■ D*£v 
decays with a small crossfeed contribution from 
Bg —)• D**£v decays. 


The continuum background is estimated using off- 
resonance data, and the B ^ Dg K£u background is 
estimated from MC simulation. We use the kinematic 
properties of the reconstructed decay to determine the 
normalisations of the other three components from data. 
For this, we consider the lepton momentum in the CM 
system of the e+e” collision, and the variable 


P*B. 


( 2 ) 


where is the energy of the Bg meson in the CM 
system approximated by \/sj2; p*g is the momentum 
of the Bg meson in the CM system approximated by 
Y^s/d — TObJ of recon¬ 

structed energies in the CM system and = \p^ +1^^ \ 
is the absolute value of the sum of the reconstructed lep¬ 
ton and Dg momenta in the CM system. When the Dg 
meson and the lepton candidate stem from the decay of 
the same Bg meson, Amis takes values larger than — 1 
because the momentum of the unreconstructed Bg decay 
products, Potheri constrained by the triangle inequality 

P*B, -P*DJ < Pother ^nd by p^^her < Kther = 

We divide the data samples into three regions: 


Amis ^ li 


B: Amis > -1 and p} < 1.4 GeV, 


C: Amis > -1 and p| > 1.4 GeV. 

As these regions are later used to determine the signal 
yields, we refer to them as “counting regions” in the 
following. Region A contains only opposite-B/^g) back¬ 
ground and can be used to determine the normalisation 
of this background. The normalisation of the other two 
components can be extracted from the measured yields 
in regions B and C, which have an enhanced fraction 
of same-B(^g) background and signal events, respectively. 
The boundary p| = 1.4 GeV is chosen to achieve ap¬ 
proximately equal event yields in regions B and C. The 
analysis is insensitive to the modeling of the Amis distri¬ 
bution for signal decays, which depends on the mass of 
the B* meson, tobj, and is thus strongly influenced by 
the poor precision on tub* ■ The semi-inclusive branching 
fractions are obtained from the relation 


BiBg dI*^X£v) 


N • 


' ^b^bX^d 


(») 


(3) 


where Agig is the measured signal yield, e is the average 
signal efficiency and is the branching fraction of 

the reconstructed Dl ' decay mode: 


Bd, = B{D- ^ 07r-; K+R -), (4) 


BDt=B{Dl- ^D:^)-Bd^. (5) 


EVENT SELECTION 

We select tracks originating from the interaction region 
by requiring \dz\ < 2.0 cm and dr < 0.5 cm, where dz and 
dr are the impact parameters along the e®" beam and in 
the transverse plane, respectively. Kaon or pion hypothe¬ 
ses are assigned to the tracks based on a likelihood com¬ 
bining the information from the Gherenkov light yield 
in the ACG, the time-of-flight information of the TOF 
and the specihc ionization dE/dx in the CDG. The kaon 
(pion) identification efficiency for tracks with a typical 
momentum of 0.75 GeV is about 96% (92%), while the 
rate of pions (kaons) being misidentihed as kaons (pions) 
is 7% (2%). The kaon and pion candidates are used to 
reconstruct Dg mesons in the high-purity decay channel 
D~ —>• (/)7r-; (j) —>• K^K~. A Dg candidate is retained 
in the analysis if it has a reconstructed mass, Mkk-k, 
within a ±65 MeV window around the nominal Dg mass, 
mo, = 1968.5MeV [31], that includes large enough side¬ 
bands to determine the combinatorial background of ran¬ 
dom KKn combinations. The reconstructed di-kaon in¬ 
variant mass, Mkk, is required to be in the mass win¬ 
dow between 1004 and 1034 MeV, corresponding to three 
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times the FWHM of the reconstructed ({> mass peak. To 
suppress combinatorial background, we impose the crite¬ 
rion I cos^heil > 0.3 on the helicity angle, defined as the 
angle between the momentum of the Dg and the K~ in 
the rest frame of the (p resonance. 

The Dg candidates with a reconstructed mass, Mkktt, 
within the range between 1953.5 and 1983.5 MeV, corre¬ 
sponding to three times the RMS of the Dg mass peak, 
are utilized for the reconstruction of D* candidates in 
the dominant decay channel —)■ Dg'y. Photon candi¬ 

dates are reconstructed from ECL clusters that are not 
attributed to a track candidate. The photon candidate 
must have a minimum energy of 125 MeV in the lab frame 
and the ratio of the energy deposit in the central 3x3 
cells of the ECL cluster to the energy deposit in the cen¬ 
tral 5x5 cells must be at least 90%. To veto photons 
from 7r° decays, we combine the photon candidate with 
any other photon candidate in the detector and require 
that the invariant mass of the two photons differs from 
the nominal 7r° mass [31] by more than 5 MeV, corre¬ 
sponding to about 0.8 times the experimental resolution. 
The angle between the Dg meson and the photon in the 
lab frame is typically less than 90°, so only candidates 
fulfilling this requirement are retained. The D* candi¬ 
dates whose mass difference between the reconstructed 
Dg and Dg candidates, AM = Mkrtt-j — Mkr^^, lies 
between 78.8 and 208.8 MeV are retained. 

Electron and muon candidates are reconstructed from 
tracks that are not used for the Dg reconstruction. Elec¬ 
trons are selected based on the position matching be¬ 
tween the track and the ECL cluster, the ratio of the 
energy measured in the ECL to the charged track mo¬ 
mentum, the transverse ECL shower shape, specific ion¬ 
ization in the CDC and the ACC light yield. Muons are 
identified using their penetration depth and the trans¬ 
verse scattering in the KLM. Hadron tracks misidenti- 
fied as leptons and leptons from secondary decays tend 
to have lower momenta than primary leptons and are 
suppressed by rejecting lepton candidates with a mo¬ 
mentum in the lab frame below 900 MeV. The electron 
(muon) identification efficiency in the selected momen¬ 
tum region is better than 89% (82%) and the proba¬ 
bility that a charged pion or kaon track is misidenti- 
fied as an electron (muon) is below 1% (2%). Leptons, 
£+, from J/tp —>• decays are vetoed by requiring 

\Mi-irh- — < 5MeV, where is the invariant 

mass of the lepton and any accepted track of the oppo¬ 
site charge, h ~, to which we assign the £ mass hypothesis. 
Furthermore, electrons are rejected if they are likely to 
stem from photon conversions, | < 100 MeV, or 

Dalitz 7r° decays, — m^oj < 32 MeV. 

We form a signal Bg candidate by pairing a can¬ 

didate with an oppositely charged lepton candidate 
To suppress background from cc continuum, we reject 
events where the normalised Dg momentum, x{Dg) = 


p* {Dg)/ ^s/4 — , is larger than 0.5 (for explanations, 

see Ref. [18]). Further suppression of the cc continuum 
background is achieved by rejecting events with a jet¬ 
like topology characterised by [cos^thrustj > 0.8, where 
fi^thrust is the thrust angle defined by the two thrust axes 
maximizing the projection of the momenta of the tracks 
and photon candidates of the 5° candidate and the rest 
of the event, respectively. 

After applying the selection criteria, 7.9% (0.4%) of the 
events contain more than one (two) Df£~ candidate(s). 
We perform a fit to the vertex of the three tracks 
used for Dg reconstruction and select the candidate with 
the best goodness-of-ht. This approach selects a correct 
candidate in 80% of the cases. The selected Dg candi¬ 
date in an event is used for £>* reconstruction; in 36.2% 
(9.7%) of the events, more than one (two) Dg'y combina¬ 
tions meet the D* requirements. We choose the photon 
candidate with the highest energy fraction deposited in 
the central 3x3 cells of a 5 x 5 cell ECL cluster. In 
the case that more than one photon candidate deposits 
all of its energy in the central 3x3 cells of the cluster, 
the candidate with the higher energy in the lab frame 
is selected. If two or more lepton candidates pass all of 
these selection criteria (2.1% of all events), we choose a 
random lepton candidate. 


FIT RESULTS 

Da fits 

We determine the yields of correctly reconstructed Dg 
mesons with binned extended maximum likelihood hts 
to the reconstructed Dg mass, M = Mrrj^, in 50 equal 
bins, indexed by j. The probability density function 
(PDF) of correctly reconstructed Dg mesons, VsigiM), 
is modeled by the sum of two Gaussian functions with 
a common mean. The PDF of the combinatorial back¬ 
ground, Vhkg{M), is a first-order Chebychev polynomial. 
We do not determine the shape parameters from simula¬ 
tion, but rather allow them to vary as free parameters in 
the fit. The Dg mass fits are performed simultaneously in 
the three counting regions {i = A,B,C) defined above. 
The width of the first Gaussian function, cti, the ratio 
of the widths of the two Gaussian functions, and the 
ratio of the normalisations of the two Gaussian functions, 
r^v, are common fit parameters in all three regions. The 
means of the Gaussian functions, pi, and the slopes of 
the polynomials describing the background, are fit¬ 
ted in each counting region individually. The likelihood 
function is: 

Tl • 

L(.-y „“»,«)= n n . (6) 

i=A.B.C ’• binsj 
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where r'c®) and 

is the vector of signal and background yields in the three 
counting regions, 9 = (cti, r^r, rjv, Ma, MS; Me, ^A, &c) 

are the shape parameters for the signal and background 
PDFs, and riij and Uij are the observed and expected 
event yields in bin j of counting region i, respectively, 


with rit = J2i ’’T'ij and i^i = 


Vi i . 

3 *-5 


The expected event 


yield, Vij, is a function of and 0: 


Mi 


-- J [vf^VsisiM) + dM. (7) 


Figures 1 (a) and (b) show the KKtt mass distributions 
together with the fit results. 


D* fits 


The D* yields are determined from binned extended 
maximum likelihood fits to the mass difference AM in 
25 equal bins, indexed by j. The combinatorial back¬ 
ground is modeled by a third-order Chebyshev polyno¬ 
mial, 'Pbkg(AM), whose parameters are constrained to 
the values obtained from fits to simulated background 
distributions. Since the background shapes vary for the 
different counting regions, the shape parameters are de¬ 
termined for each counting region separately. The signal 
peak is modeled by the sum of a Gaussian function and a 
Crystal Ball function [42] to account for energy loss due 
to material in front of the calorimeter: 


7^sig(AM) oc 


fexp(-l^) 
+ \ 

I ^ AM-,j. y 


if > -a 

(T 

■ c AM-u , 
if < -a 


A common mean, fi, is used for both the Gaussian and 
the Crystal Ball functions. We perform a fit to the sim¬ 
ulated signal distribution and fix the parameters tat, To-, 
a and n at the obtained values. The width a and the 
mean of the signal peak /i are varied in the fit to data; 
the parameter tr is fitted simultaneously in all counting 
regions while \i is fitted individually for each counting 
region. The likelihood function is constructed analogous 
to Eqs. (6) and (7) with additional factors, to implement 
the constraints of the background PDF parameters tak¬ 
ing into account their correlations. The results of the 
AM fits in the different counting regions are presented 
in Figs. 1 (c) and (d). 


the off-resonance data sample is not sufficient to deter¬ 
mine the shape parameters in the fits, they are fixed to 
the values obtained in the fits to T{5S) data in the cor¬ 
responding counting region. The CM energy, y/s, in the 
expression for the Amis variable in Eq. (2) is replaced by a 
constant value of 10.876 GeV because, otherwise, the de¬ 
nominator would not be defined. The cc continuum yields 
from the fits to off-resonance data are multiplied by the 
scale factor S = {Cr{ 5 S)/sr{ 5 S))/{^os/sos) = 1.81±0.02 
to account for the differences in integrated luminosities, 
C, and the 1/s dependence of the e+e“ —)■ cc cross sec¬ 
tion. Additionally, a shape correction for differences of 
the yields in the counting regions between off-resonance 
and T(5S') data is determined from MC simulation and 
applied. The small background from B ^ Dl 'K£v de¬ 
cays is estimated from MC simulation using a simple 
phase space model. The backgrounds from continuum 
processes and B ^ Dl ' K£v decays are subtracted in 
each counting region from the yields measured in T(5S') 
data. 


Signal extraction 

After subtraction of the continuum and the B —>■ 
Di K£v background components, the remaining yields 
contain three contributions: opposite-B(^s^ background, 
same-B(^g) backgrounds and signal. The three contribu¬ 
tions are constrained by the event yields in the three 
counting regions. We introduce a scale factor, aj , for each 
contribution, j. The determination of the scale factors is 
equivalent to solving a system of three linear equations 
with three unknowns. In order to obtain the uncertain¬ 
ties on the scale factors, we minimize: 

= V ~ 

where the index i runs over the three counting regions, 
Di is the event yield determined by the fits to the Mkkh 
or AM distributions in data, and AD^ is the statistical 
uncertainty of these fits, Ni j is the MC prediction for the 
contribution j, and is its statistical uncertainty. 

Table II lists the scale factors, aj, obtained from the 
minimization and the signal yields. 


Background subtraction 






(9) 


To estimate the cc continuum background, the Dg and 
D* yields are measured in and samples re¬ 

constructed in the off-resonance data. Since the size of 


Figure 2 shows the Amis and p*{t) distributions in the 
three counting regions after applying the scale factors 











Events / 5.2 MeV Events / 5.2 MeV Events / 2.6 MeV Events / 2.6 MeV 
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{a.)Ds e+ 





(h)D7li+ 






{c)Dt-e+ 





{d)Dt-ii+ 


FIG. 1: The Mkktt distributions for D~£'^ events and AM distributions for D*~£~^ events reconstructed in the T(5S) data 
for the three counting regions. The black points with uncertainty bars are the data, the red solid curve represents the total fit 
result, and the green dashed line is the fitted background component. 



























TABLE II: The scale factors, aj, for the MC components obtained by minimizing the function defined in Eq. (8). The errors 
are the statistical uncertainties of the data and the MC sample. The signal yields are determined from Eq. (9). The yields of 
the other components are given in Tables III and IV. The signal efficiencies are obtained by averaging over the efficiencies for 
the Dstv, D*£u and D**lv modes, taking into acconnt the expected relative abundance of the signal components. The given 
errors of the efficiencies are the statistical nncertainties of the MC sample. 


Channel 

Scale factors 

Opposite-B^g) Same-Bf^g-) Signal 

Signal yield Efficiency [%] 

DgXeu 

DgXpu 

D*gXev 

D*gXp.v 

1.02 ± 0.04 1.00 ± 0.20 1.06 ± 0.04 

1.06 ± 0.04 0.94 ± 0.16 1.09 ± 0.04 

0.89 ± 0.12 1.66 ± 0.71 1.00 ± 0.11 

0.96 ± 0.12 1.50 ± 0.58 1.13 ± 0.12 

4470 ± 161 16.9 ± 0.1 
4411 ± 161 16.3 ± 0.1 

724 ± 79 4.6 ± 0.1 

804 ± 86 4.6 ± 0.1 



FIG. 2: Distributions of Xmis and p*{£) for reconstructed oi*'' C*" events. The black points with uncertainty bars show the 
yields in the T(5S) data determined by fits to the Mkk-w distributions for D~ and the AM distributions for . 

The stacked histograms represent the signal and background expectations after applying the scale factors aj (see Table II). The 
components are, from bottom to top: continuum background (white), B —>■ Oi*-’Klu background (dark green), opposite-B^s) 
primary leptons (solid blue), opposite-B^g'j secondary leptons and misidentified hadrons (hatched blue), same-Bf^gf background 
(hatched yellow), signal Bg —>■ Dg£u (solid red), signal Bg —>■ Dgiv (hatched red), signal Bg —>■ D**lv (cross-hatched red). The 
vertical black line illustrates the division of the counting regions. The displayed binning of the Xmis and p* (i) distributions is 
used only to illustrate the data-MC agreement; the signal yield, Asig, is extracted from the measured Di*^~ yields in the three 
counting regions A, B and C listed in Tables III and IV. 
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TABLE III: The Dj£'^ yields obtained from the Mkktt fits to T(5S') data in the three counting regions (A, B, C) and the 
corresponding signal and background expectations. The scale factors from Table II obtained by minimizing Eq. (8) are applied 
to the MC expectations (3) - (5). The errors are the statistical uncertainties of the data and MC samples, respectively, and do 
not contain the scale factor uncertainties. Uncertainties are omitted if they are smaller than 0.5. 




Electrons 



Muons 



A 

B 

C 

A 

B 

C 

T(5S') data 

1807 ± 53 

4274 ± 87 4215 ± 82 

1902 ± 54 4544 ± 89 4375 ± 81 

(1) Continuum (scaled off-resonance data) 

130 ± 34 

278 ± 37 

137 ± 22 

102 ± 32 

298 ± 40 

134 ± 25 

(2) B DsK£v 

0 

48 ± 7 

18 ± 4 

0 

46 ± 7 

18 ± 4 

(3) Opposite-B( 3 ), secondary leptons, mis-ID hadrons 

110 ± 4 

555 ± 10 

61 ± 3 

205 ± 6 

826 ± 12 

107 ± 4 

(3) Opposite-i3(3), primary leptons 

1565 ± 16 

1165 ± 14 

1032 ± 13 

1594 ± 17 

1081 ± 14 

1043 ± 14 

(4) Same-B(s) background 

0 

638 ± 10 

89 ± 4 

1 

798 ± 11 

158 ± 5 

(5) Signal {Dglv) 

0 

492 ± 9 

669 ± 11 

0 

489 ± 9 

693 ± 11 

(5) Signal {D*lv) 

1 

951 ± 13 

2072 ± 19 

0 

872 ± 13 

2072 ± 19 

(5) Signal {DT(.u- DT ^ Dl) 

0 

28 ± 2 

41 ± 3 

0 

26 ± 2 

40 ± 3 

(5) Signal {DTiv, DT ^ Dt) 

0 

117 ± 5 

98 ± 4 

0 

109 ± 4 

110 ± 4 


TABLE IV: The D*~£'^ yields obtained from the AM fits to T(55) data in the three counting regions (A, B, C) and the 
corresponding signal and background expectations. The scale factors from Table II obtained by minimizing Eq. (8) are applied 
to the MC expectations (3) - (5). The errors are the statistical uncertainties of the data and MC samples, respectively, and do 
not contain the scale factor uncertainties. 




Electrons 



Muons 



A 

B 

C 

A 

B 

C 

T(5S') data 

336 ± 33 

656 ± 48 

662 ± 46 

370 ± 35 

739 ± 52 

741 ± 50 

(1) Scaled off-resonance data 

32 ± 22 

61 ± 17 

24 ± 11 

49 ± 19 

54 ± 18 

20 ± 11 

(2) B DsKtv 

0 

6 ± 2 

2 ± 1 

0 

4 ± 2 

2 ± 1 

(3) Opposite-R(s), secondary leptons, mis-ID hadrons 

24 ± 2 

60 ± 3 

4 ± 1 

48 ± 3 

99 ± 4 

13 ± 1 

(3) Opposite-R(s), primary leptons 

279 ± 6 

147 ± 5 

120 ± 4 

273 ± 7 

147 ± 5 

109 ± 4 

(4) Same-R( 5 ) background 

0 

151 ± 6 

20 ± 2 

0 

188 ± 7 

39 ± 3 

(5) Signal {Dllv) 

0 

227 ± 6 

483 ± 9 

0 

241 ± 7 

547 ± 10 

(5) Signal {DT£v, DT ^ Dl) 

0 

6 ± 1 

8 ± 1 

0 

6 ± 1 

11 ± 1 


SYSTEMATIC UNCERTAINTIES 

The different sources of systematic uncertainties on the 
measured signal yields are described below. They com¬ 
prise detector effects and the modeling of the signal and 
backgrounds. An overview can be found in Table V. 

Detector effects 

The uncertainty on the track finding efficiency is 0.35% 
per track and thus 1.4% for four tracks. The pho¬ 
ton efficiency is studied with radiative Bhabha events, 
from which the uncertainty is estimated to be 2%. The 
calibration of kaon and pion identification efficiencies 
is estimated from a sample of reconstructed D*'^ —>■ 
ZI°7r+; —>■ K~t:^ decays. A variation of the obtained 
calibration factors within their uncertainties changes the 
measured signal yield by 1.4%. The efficiency of the lep¬ 
ton identification is estimated using the two processes 
77 —>■ £+£“ and J/ijj —t ■ The corresponding uncer¬ 

tainties on the measured signal yields are 1.0% and 1.6% 
for the electron and muon modes, respectively. The rates 
of hadrons being misidentified as leptons are estimated 


with the aforementioned D*^ sample. The uncertainties 
due to this estimation are 0.1% {DgXeiy), 1.3% {D^X^v), 
0.1% {DlXev) and 1.9% {DlX^iv). 


Signal and background modeling 

To study uncertainties of the PDFs in the and D* 
fits, we repeat the fits with alternative fit models and as¬ 
sign the resulting change of the signal yield as systematic 
uncertainty. Herein, we focus on the tails of the signal 
peaks because they can be easily assigned in the fit to the 
background component without deteriorating the agree¬ 
ment of the data with the fitted curve. The signal PDF 
in the Dg fits is modified by replacing the second Gaus¬ 
sian function by a bifurcated Gaussian function. This 
choice is motivated by a small asymmetry of the signal 
peak due to final state radiation. The normalization and 
the widths of the bifurcated Gaussian function are de¬ 
termined relative to the normalization and width of the 
Gaussian function from a fit to signal MG. These param¬ 
eters are fixed in the fit to data. Based on the observed 
change of the signal yield, we assign a 3% PDF uncer¬ 
tainty. In the D* fits, the tails of the signal peak are 
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TABLE V: Relative systematic uncertainties on the signal 
yields in %. 


DgXev DgXgi/ D*Xev D*Xp,i/ 

Detector 





Tracking efficiency 

1.4 

1.4 

1.4 

1.4 

Photon efficiency 

— 

— 

2.0 

2.0 

Kaon and pion ID 

1.4 

1.4 

1.4 

1.4 

Lepton efficiency 

1.0 

1.6 

1.0 

1.6 

Hadron misidentification 

0.1 

1.3 

0.1 

1.9 

Signal and background modeling 




PDF for Mkktt and AM fits 

3.0 

3.0 

5.0 

5.0 

Continuum shape 

1.2 

0.3 

1.2 

0.3 

B Kin modeling 

0.3 

0.3 

0.1 

0.1 

Signal 





Composition 

4.8 

4.8 

0.3 

< 0.1 

Form factors 

0.9 

1.0 

1.0 

1.0 

Efficiency 

3.1 

3.1 

3.0 

3.0 

Opposite-B(^s'j background 





Composition 

1.6 

2.2 

1.0 

2.5 

Bg fraction 

0.2 

0.2 

< 0.1 

< 0.1 

Shape 

1.0 

1.0 

1.0 

1.0 

Same-B(g) background 





Composition and shape 

0.3 

0.3 

0.4 

0.7 

Bg production mode 

0.1 

0.1 

0.3 

0.3 

Beam energy 

1.0 

1.0 

0.5 

0.5 

Total 

7.3 

7.6 

6.9 

7.6 


described by the Gaussian component of the signal PDF. 
When this Gaussian function is removed from the sig¬ 
nal PDF, i.e. a Grystal Ball function only is used (cf. 
Ref. [43]), the signal yields decrease by 5%. Hence, we 
estimate the PDF uncertainty with 5%. 

The uncertainty due to the continuum scale factor, S, 
is negligible. The uncertainty due to the shape correc¬ 
tion for the continuum background is estimated as the 
full difference of the result with and without the correc¬ 
tion applied, which is 1.2% and 0.3% for electrons and 
muons, respectively. To estimate the influence of the 
choice of the B ^ Dg ’ Kiv decay model, we replace the 
phase space model used in the nominal result with the 
ISGW2 model [44], assuming that the decay proceeds via 
B —)■ Dq —> Di*'^K. The use of this alternative 

model increases the signal yields by 0.3% and 0.1% for 
the DgXlv and DlXt-v channels, respectively. We also 
vary the B ^ Dl Klv branching fraction by the mea¬ 
sured uncertainty and observe no significant change in 
the measured yield. We test the stability of the signal 
extraction when the boundary between counting region 
B and C is varied between p*^ = 1.3 and 1.5 GeV. The 
resulting change of the signal yields is consistent with the 
expected change due to the increase/decrease of statis¬ 
tics in the respective counting regions and, therefore, no 
systematic uncertainty is assigned. 

The systematic uncertainty on the signal composition 
in the ngXiv channels is obtained by evaluating the ef¬ 
fect of scaling the relative amount of Bg —> Dlii/ de¬ 


cays up and down by 30% and adjusting the Bg —>■ Dgiv 
component such that the total number of MC events is 
conserved. This variation covers most of the recent the¬ 
ory predictions and causes a 4.7% change of the signal 
yields. To estimate the impact of D** —>• DgX crossfeed, 
we double the Bg —>■ 0 ** 11 / contribution in the signal 
component, which increases the signal yield by 1%. The 
D^Xiv signal component is expected to be dominated 
by Bg —>• D*£i/ decays and hence the uncertainty due to 
the amount of D^* D*X crossfeed is negligible for this 
channel. 

The B — D^*'>ii/ form factor parameters from Ref. [37] 
used to simulate the Bg —>■ Dg ' iv decays are measured 
with an accuracy of 2-3%. However, SU(3) flavor sym¬ 
metry breaking effects may cause deviations at the order 
of 10% [4]. To account for these differences, we vary 
each form factor parameter of a given decay indepen¬ 
dently up and down by 10%. The resulting average de¬ 
viation from the nominal signal yield is added linearly 
for each variation. The uncertainty of the LLSW model 
for Bg —>• D**£i/ decays is evaluated by repeating the 
measurement with different sets of model parameters, as 
specified in Ref. [38]. The total systematic uncertainty 
due to form factor modeling is given by the quadratic 
sum of the uncertainties from all decay modes and does 
not exceed 1%. 

The signal efficiencies are studied in bins of three dis¬ 
tributions: the lepton momentum, the Dg momentum, 
and the angle between the reconstructed Dg meson and 
the lepton in the CM system. A re-calculation of the 
average efficiencies based on the observed data yields 
changes the signal by at most 3.1%. 

The modeling of the opposzte-H(s) component is stud¬ 
ied in same-sign control samples. The same- 

sign selection ensures that these samples contain only 
opposite-Bf^g) combinations. Compared to the D\ t*' 
samples, the relative contribution of Bg decays is en¬ 
hanced in this control sample. Two components of 
the opposite-B(^g'j sample are distinguished: (1) primary 
leptons and (2) secondary leptons and hadron tracks 
misidentified as leptons. Scale factors for the normali¬ 
sation of these two MC components are determined from 
fits to the p*{i) distributions of the samples. The 

obtained scale factors are in agreement within the fit un¬ 
certainties of about 10%. A variation of the normali- 
sations of the two components in the Dg samples 
within this 10% uncertainty changes the signal yields be¬ 
tween 1.0% and 2.5%, depending on the reconstructed 
channel. We also vary the fraction of Bg decays in the 
opposite-B(^g'f component by 20%, corresponding to the 
uncertainty of the Bg production rate, fg [31]. The re¬ 
sulting change of the signal yields is less than 0.2%. The 
shape uncertainty of the oppos*te-i?(s) component is eval¬ 
uated in a data-driven way by using again the 
samples, from which the event yields are determined in 









11 


the three counting regions with the identical procedure 
as applied in the measurement. We calculate the ratios 
of data and MC yields for each counting region. These 
ratios range from 0.86 to 0.91 for electrons and from 0.96 
to 0.97 for muons. We then modify the MC predictions 
for the opposite-Bf^g^ component in the Dji~^ simulation 
accordingly and study the impact on the measurement. 
The results change by less than 0.4%, so an uncertainty 
of 1% on the modeling of the opposite-B(^g-^ component 
is a reasonable estimate, considering the differences be¬ 
tween the control samples and the D~t^ signal 

samples. The described approach cannot be transferred 
to the D*~t^ measurements because of the smaller sam¬ 
ple sizes. However, the composition of the opposite-B(^g) 
background in the D*g~t^ sample is similar to the one 
in the Dj sample and hence the same uncertainty is 
assigned. 

The decays contributing to the same-B(^g'^ background 
component can be grouped into four classes with the cor¬ 
responding fraction in the electron/muon channel given 
in parentheses: i?(s) —decays (70% / 48%), lep¬ 
tons stemming from r produced via Bg and Dg decays 
(21% / 16%) and hadrons misidentified as leptons (9% / 
34%). There are no significant differences in the compo¬ 
sition between the Dj and the D*~£~^ channels. We 
vary the fraction of leptons from r decays and the frac¬ 
tion of misidentified hadrons by ±50% and take half the 
difference of the resulting signal yields as the systematic 
uncertainty, which is below 1% for all measurements. Po¬ 
tential modeling uncertainties of the same-B(^s^ compo¬ 
nent are assumed to be covered by the large variation of 
the composition. 

We estimate the impact of the uncertainty on the dif¬ 
ferent Bg production channels at the T(5S') energy by 
scaling the BgB* component up and down by 30% and 
assign half of the change in the signal yield as the sys¬ 
tematic uncertainty of 0.1% and 0.3% for the D~i'^ and 
Dg~£~^ modes, respectively. The beam energy is con¬ 
servatively varied by ±3 MeV and signal yield variations 
of 1% and 0.5% are observed for the Dj£'^ and 
modes, respectively. 


RESULTS AND DISCUSSION 

The semi-inclusive semileptonic Bg branching fractions 
are calculated from Eq. (3). Since the Dg 
reconstruction mode is also used in the determination of 
Nb^Bs [2^] I ttie B{Df —>• Ar+A“7r+) branching fraction 
cancels out. Using the branching fraction ratio B{D'^ —> 
(j)TT^)/B{Dg —)• A'+A“7r+) = (41.6 ± 0.8)% and the 


branching fraction B{D*g —)• Ugy) = (94.2 ± 0.7)% [31], 
we obtain the semi-inclusive branching fractions: 

DgXev : [8.1 ± 0.3(stat) ± 0.6(syst) ± 1.4(ext)]% , 

DgXpLv : [8.3 ± 0.3(stat) ± 0.6(syst) ± 1.5(ext)]% , 

D*gXev : [5.2 ± 0.6(stat) ± 0.4(syst) ± 0.9(ext)]% , 

D*gXpu : [5.7±0.6(stat)±0.4(syst)±1.0(ext)]%. 

The first uncertainty is the statistical uncertainty of the 
data and MC samples, the second is the systematic un¬ 
certainty of the measurement, and the last uncertainty is 
due to the external measurements of Nb^b^ Bj^m. 
The electron and muon samples are statistically indepen¬ 
dent because only one candidate is selected per event. 
Taking into account that the systematic uncertainties 
are all correlated except the one for lepton identifica¬ 
tion, we calculate the combination of the measurements 
as weighted averages: 

DgX£v : [8.2 ± 0.2(stat) ± 0.6(syst) ± 1.4(ext)]% , 
D*gX£v : [5.4±0.4(stat)±0.4(syst)±0.9(ext)]%. 

The obtained Bg DgX£v branching fraction can be 
compared to the difference between the inclusive branch¬ 
ing fraction, B{Bg —)■ X(.£v), and the branching frac¬ 
tion of the D*g*£v modes, where the D** does not de¬ 
cay to a Dg meson. The value of B{Bg —>■ Xc£v) is es¬ 
timated to be (10.0 ± 0.4)%, using the branching frac¬ 
tion B{B^ —>■ Xc£v) [31, 46], an estimate for the ra¬ 
tio of the semileptonic widths of the Bg and B^ meson, 
rsi(7?s)/rsi(H°) = 0.99 [47] and the measured lifetimes 
of the and Bg mesons [31]. 

We assume that only the semileptonic decay modes 
with Zlsi(2536) and Ils2(2573) mesons do not contain 
Dg mesons in the final state. We obtain the estimate, 

BestiBg ^ DgX£u) 

= B{Bg ^ XM ■ [1- 

B{Bg Dg 2 X£iy)/B{Bg XJv)- (10) 

B{Bg ^ Dg^X£v)IB(Bg ^ Xdv)] 

= (9.1 ±0.4)%, 

where the ratios B{Bs —>■ Dg 2 X£v) /B{Bg 

Xciv) = [3.3 ± l.O(stat) ± 0.4(syst)]% and B{Bg —> 
DgiX£v)IB{Bg XJv) = [5.4±1.2(stat)±0.5(syst)]% 
were measured at LHCb [45]. The result of our mea¬ 
surement is in agreement with the estimate, BestiBg —>■ 
DgX£v). The rate of Bg D*g*£v-D*g* D*gX 

decays can be constrained from the comparison be¬ 
tween the measured semi-inclusive branching fraction, 
B{Bg —>■ D*gX£v) with the exclusive theory predictions 
for B{Bg —)• D*g£v). For example, using, the predic¬ 
tion from Ref.[8], one obtains B{Bg D**£iy;Dg* —> 
DgX) < 2.0% at the 90% confidence level. 
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The measurement can also be used to determine Ng g 
using the estimate of the Bs —>■ DgXiv branching frac¬ 
tion from Eq. (10): 


N,iJ[eB{Ds ^ <j>{K+K-)n+)] 
2Best(Bs ^ DsXev) 


( 11 ) 


For Nsig/e, we insert the weighted average of the elec¬ 
tron and muon modes, -/Vgig/e = [26.7 ± 0.7(stat) ± 
2.0(syst)] X 10^; for B{Ds —)• we use 

the value (2.24 ± 0.10)% [31]. We obtain = 

[6.53 ± 0.17(stat) ± 0.49(syst) ± 0.41(ext)] x 10®, corre¬ 
sponding to the cross section (T(e+e“ —>■ = 

[53.8 ± 1.4(stat) ± 4.0(syst) ± 3.4(ext)] pb at the CM en¬ 
ergy -y/s = 10.86 GeV. The first two uncertainties are 
the statistical and systematic uncertainties from the mea¬ 
surement, respectively, and the last uncertainty is due to 
B{D, (j){K+K-)Tr+) and BestiBs ^ DsXlv). The ob¬ 
tained result is in agreement with Ng g = (7.1 ± 1.3) x 
10® obtained by Belle with a different technique [23] and 
has a significantly improved precision. 


SUMMARY 

We have presented the first measurements of the 
semi-inclusive branching fractions of Bs —>■ DsXii/ and 
Bg —)■ D*Xii' decays. The measured branching fractions 
are B{Bs —>■ DgXiv) = [8.2 ± 0.2(stat) ± 0.6(syst) ± 
1.4(ext)]% and B{Bs D^Xiv) = [5.4 ± 0.4(stat) ± 
0.4(syst) ± 0.9(ext)]%. In addition, the analysis of these 
decays provides the currently most precise estimate of 
the Bl 'Bs production cross-section at the CM en¬ 
ergy = 10.86GeV: cr(e+e- ^ Bi*'’Bi*^) = (53.8 ± 
1.4(stat) ± 4.0(syst) ± 3.4(ext))pb. 
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